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ABSTRACT 



Based on early solar system abundances of short-lived radionuclides (SRs), 
such as ^'^Al (Ti/2 = 0.74 Myr) and ^°Fe (T1/2 = 1.5 Myr), it is often asserted that 
the Sun was born in a large stellar cluster, where a massive star contaminated the 
protoplanetary disk with freshly nucleosynthesized isotopes from its supernova 
(SN) explosion. To account for the inferred initial solar system abundances of 
short-lived radionuclides, this supernova had to be close (~ 0.3 pc) to the young 

1 Myr) protoplanetary disk. 

Here we show that massive star evolution timescales are too long, compared 
to typical timescales of star formation in embedded clusters, for them to explode 
as supernovae within the lifetimes of nearby disks. This is especially true in an 
Orion Nebular Cluster (ONC)-type of setting, where the most massive star will 
explode as a supernova ~ 5 Myr after the onset of star formation, when nearby 
disks will have already suffered substantial photoevaporation and/or formed large 
planetesimals. 

We quantify the probability for any protoplanetary disk to receive SRs from a 
nearby supernova at the level observed in the early solar system. Key constraints 
on our estimate are: (1) SRs have to be injected into a newly formed 1 Myr) 
disk, (2) the disk has to survive UV photoevaporation, and (3) the protoplanetary 
disk must be situated in an enrichment zone permitting SR injection at the solar 
system level without disk disruption. The probability of protoplanetary disk 
contamination by a supernova ejecta is, in the most favorable case, 3 x 10^"^. 



Subject headings: solar system: Origins, Sun: Short-lived radioactivities; Super- 
novae; Meteorites; ISM: individual (Carina Nebula), stars: formation 
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1. Introduction 

Short-lived radionuclides (SRs) are radioactive elements with half-lives of the order of 1 
Myr. Their presence in the nascent solar system is inferred from an excess of their daughter 
isotopes in meteorite components (Russell, Gounelle & Hutchison 2001). Some SRs [^'^Be 
(Ti/2 = 1.5 Myr), 26A1 (T1/2 = 0.74 Myr), ^^Cl (T1/2 = 0.30 Myr), ^iQa (T1/2 = 0.10 Myr), 
^^Mn (Ti/2 = 3.7 Myr), ™Fe (T1/2 = 1.5 Myr)] were present in the protoplanetary disk at 
abundances substantially higher than the levels expected for average interstellar medium 
(Meyer & Clayton 2000). These SRs cannot therefore be inherited from the interstellar 
medium and require a last minute origin (e.g. Wadhwa et al. 2007). The origin of SRs is 
highly debated in the cosmochemistry community, because it has important consequences 
for e.g. early solar system chronology (e.g. McKeegan & Davis 2004; Gounelle & Russell 
2005a,b) and planetesimal heating (Urey 1955; Ghosh & McSeen 1998). In addition, under- 
standing the origin of SRs will provide key constraints on the astrophysical environment in 
which our solar system was born (e.g. Hester & Desch 2005; Goswami et al. 2005; Gounelle 
2006). 

Generally speaking, SRs can be made either by thermal nucleosynthesis, at KeV energies 
in the interior of stars, or by non-thermal nucleosynthesis, i.e. nuclear reactions involving 
cosmic rays or accelerated particles at MeV energy from the Sun. There are therefore two 
possible sources for short-lived radionuclides in the solar protoplanetary disk: (a) Injection 
into the molecular cloud core or protoplanetary disk from a nearby, late-type star (e.g. Busso, 
Gallino & Wasserburg 2003), (b) In situ production in the protoplanetary disk by irradiation 
of dust and gas with accelerated particles from the proto-Sun (e.g. Gounelle et al. 2006). 

We focus here on the possible delivery of SRs by a late type star. Asymptotic Giant 
Branch (AGB) stars are possible candidates (Wasserburg et al. 2006), but this source of SRs 
in the early solar system is considered highly unlikely because of the low probability of an 
encounter between an AGB star and a star forming region (Kastner & Myers 1994). Type II 
supernovae represent an alternative to AGB stars as a source of SRs (Cameron et al. 1995). 

Of models based on a supernova (SN) origin for SRs there are two types. Either a SN 
injects freshly synthesized SRs into a nearby molecular cloud core, triggering its gravitational 
collapse (Cameron & Truran 1977; Cameron et al. 1995; Boss & Vanhala 2000), or directly 
into a nearby protoplanetary disk (Chevalier 2000; Hester & Desch 2005; Ouellette et al. 
2005). The first SN scenario is now considered less likely because only very specific conditions 
allow a supernova Shockwave to trigger the collapse of a molecular cloud core and, at the 
same time, inject SRs (e.g. Boss & Vanhala 2000). In the second scenario, which is currently 
receiving a lot of attention (Chevaher 2000; Hester & Desch 2005; Ouellette et al. 2005; 
Ouellette, Desch & Hester 2007), the SN has to be very close (~ 0.3 pc) to the protoplanetary 
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disk in order to allow the disk to intercept enough SN ejecta to account for the solar system 
inventory of SRs. It is thus assumed that the massive star, which evolved into a SN, and 
the protoplanetary disk were coeval and formed in the same stellar cluster (e.g. Hester & 
Desch 2005). Massive stars form in large stellar clusters (Lada & Lada 2003). The nearby 
SN injection model therefore implies that our Sun was born in a large stellar cluster. The 
Orion Nebula Cluster (ONC) is often cited as a good analog for such an environment (Hester 
& Desch 2005). 

However, if the Sun was born in a stellar cluster massive enough to contain a SN, it 
has important implications for the childhood of our solar system. Massive stars arc known 
to power winds carving large bubbles of hot gas in surrounding cold molecular gas (Weaver 
et al. 1977). They also emit large UV fluxes, which create HE bhsters in cold molecular 
gas (O'Dell 2001) and photoevaporate nearby protoplanetary disks (e.g. Johnstone, Bally 
& HoUenbach 1998). The presence of numerous massive stars close to planetary systems 
can also dynamically modify the orbits of the giant planets and of Kuiper Belt Objects (e.g. 
Adams & Laughlin 2001). It is therefore essential to evaluate the a priori probability of our 
Sun to be born in such an harsh astrophysical setting. 

Adams & Laughlin (2001) estimated the probability of the Sun to have been born in a 
stellar environment large enough to contain a massive star (which could deliver SRs via a 
SN explosion) and small enough to preserve both the orbits of the giant planets as well as 
those of the Kuiper belt objects. These authors made the assumption that any star more 
massive than 25 Mq was able to deliver SRs to the solar protoplanetary disk, and they did 
not place constraints on the time window during which SRs were injected into the disk. The 
spatial structure of the cluster was not included in this model. Subsequently, WiUiams & 
Gaidos (2007) provided an estimate for the probabihty of the nearby SN injection scenario. 
However this estimate did not take quantitatively into account the disk photoevaporation and 
overestimated the time window during which SN injection of SRs in a nearby protoplanetary 
disk can occur. 

The present work is organized as follows. First, we review the cosmo chemical constraints 
on the origin of SRs (§2). Second, we present recent observations and models of molecular 
cloud dynamics, star formation rate, accretion disk lifetime and evolutionary timcscalcs of 
massive stars pertinent to an evaluation of the proposed ONC-like setting for the formation 
of the solar system (§3). Third, we calculate the probability for a protoplanetary disk 
attached to a low-mass star to receive SRs from a nearby SN at the levels inferred for the 
early solar system, with the necessary constraints that SR delivery happens within the time 
window specifled by the cosmochemical data, and that the disk is not destroyed by UV 
photoevaporation (§4). Our model assumptions and limitations are discussed in section 
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§5. In section §6 we present the implications of our findings for the origin of short-hved 
radionuchdes. 



2. Polluting our protoplanetary disk with a SN ejecta: Timing and distance 

Among the SRs that could have been delivered by a supernova to the protoplanetary 
disk around the young Sun, ^^Al provides the strongest constraint on the timing of such 
an event (e.g. Wadhwa et al. 2007; Russell et al. 2006). Aluminium-26 is more abundant 
in primitive chondrite components, such as Calcium-, Aluminium-rich Inclusions (CAIs) 
and chondrules (e.g. Galy et al. 2000) than in more evolved asteroidal crustal rocks (e.g. 
Srinivasan et al. 1999). CAIs have the oldest measured Pb-Pb ages of any solid matter in 
the solar system (Bouvier et al. 2007). Some iron meteorites, believed to sample the cores 
of differentiated asteroids, have model Hf-W ages as old as those of CAIs (Markowski et 
al. 2006). They must therefore have formed contemporaneously with CAIs. This implies 
that ^^Al, which is beheved to have been the main heat source responsible for asteroidal 
differentiation, was present in the protoplanetary disk already at the time of iron meteorite 
parent-asteroid accretion (Bottke et al. 2006). Both observations suggest that ^^Al was 
present in the protoplanetary disk from the early stages of its evolution, certainly as far 
back in time as the meteorite record can bring us (e.g. Wadhwa et al. 2007; Russell et al. 
2006). Although the exact timescale is not known, it is reasonable to assume that if ^^Al 
was delivered to the solar system from a SN source, this injection occurred during the class 
or class I stage of our protoSun (e.g. Montmerle et al. 2006). The combined duration of 
the class and class I phases is ~ 10^ yr. Here we will conservatively assume that the SN 
delivery of SRs to the protoplanetary disk occurred within a time interval, ©, of ~ 1 Myr 
after the formation of the disk (see also Ouellette et al. 2005; Looney, Tobin & Fields 2006). 

The distance r between the SN and the protoplanetary disk cannot be too small, oth- 
erwise the disk will be destroyed by the impact of the ejecta. Nor can this distance be too 
large because this would limit the amount of injected SRs. There is a restricted distance 
interval, the 'minimum' (rmin) and the 'maximum' (rmax) acceptable distance between the 
SN and protoplanetary disk, referred to as the 'enrichment zone', which satisfies these two 
important constraints. 

Chevalier (2000) estimated that at distances smaller than 0.25 pc, momentum transfer 
between the SN ejecta and the disk will cause the later to be stripped away. Based on 2D 
numerical simulations, Ouellette, Desch & Hester (2007) concluded that a disk as close as 
0.1 pc from a 25 Mq SN can survive disruption. For the minimum distance, rmin, we will 
adopt a conservative value of 0.1 pc. 
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The maximum distance rmax, below which the protoplanetary disk receives enough SRs 
to account for solar system abundances depends on a limited number of parameters. Though 
it is beyond the scope of the present paper to explicitly derive Tmax, we indicate how it can 
be calculated and summarize the estimates given by several groups (Ouellette et al. 2005; 
Looney, Tobin & Fields 2006). 

For any radionuclide, the mixing fraction, /, of the supernova ejecta mixed into the 
protoplanetary disk is: 

/ = ^ X e''/', (1) 

where Mss is the initial mass of the radionuclide in the solar system (in M©), Ysn is the 
total yield of the radionuclide from the considered supernova (in M©), r is the mean life, t] 
the injection efficiency and A is the time interval between the end of nucleosynthesis in the 
supernova and incorporation into refractory meteorite components (CAIs). Note that / <^ 1 
as Mss ^ Ysn and A ^ r (see below). Following Cameron et al. (1995), we translate the 
mixing fraction / into a geometrical relationship, r = ya/1//, where tq is the size of the 
disk and r the distance between the disk and the massive star. Putting these two expressions 
together, we obtain the general expression: 



Mss is estimated from the disk mass, the solar abundances and the initial abundance 
of the considered radionuclide identified to that of CAIs (e.g. Gounelle & Meibom 2007). 
Using SN yields calculated by state-of-the art nucleosynthetic models corresponding to a 
large range of progenitor masses and metallicities (Rauscher et al. 2002; Woosley & Weaver 
1995), and assuming an injection efficiency of 1, Looney, Tobin & Fields (2006) as well as 
Ouellette et al. (2005) calculate that the solar protoplanetary disk had to be at a distance 
between 0.02 and 0.3 pc from the massive star for disk sizes varying between 30 and 100 
AU. A is constrained by the decay of short-lived radionuclides and is of the order of ~ 1 
Myr. Note that the disk size cannot be significantly larger than 100 AU because, before any 
injection occurs, the disk has been sitting nearby a strongly ionizing O star for some Myr (see 
section [3j), resulting in severe disk truncation (Ouellette et al. 2005). The alternative case of 
larger disks examined by Looney, Tobin & Fields (2006) is therefore not relevant here. For 
more realistic values of t] (i.e. r] <^ 1) the injection distance is reduced to values approaching 
the minimum acceptable distance (0.1 pc), which would cause a problem with regard to 
disk survivability. We will consider here that a generous maximum distance between the 
protoplanetary disk and the supernova is rmax = 0.4 pc. 

Therefore, the acceptable distance interval for which the SN explosion delivers SRs at 




(2) 
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the level observed in the solar system, but does not destroy the protoplanetary disk is quite 
narrow. The enrichment zone lies therefore between r^^^ = 0.1 pc and r^ax = 0.4 pc. 

3. The unlikely ONC-like setting 

The Orion Nebular Cluster is the epitome of a massive star forming region in the solar 
neighborhood, often invoked as an analog for the birthplace of our Sun (Hester & Desch 
2005; Ouellette et al. 2005; Bally, Moeckel & Throop 2005). The ONC is dominated by 
the four bright Trapezium O and B stars, which are named 6^ and numbered from A to D 
(O'Dell 2001). It is situated ~ 450 pc away from the Sun and contains roughly half a dozen 
stars massive enough to end their lives as SN (Hillenbrand 1997). In the ONC, the most 
massive object is the 07 star 6^C Ori, which has a mass of ~ 40 Mq. Within 1 pc of 6^C 
Ori are thousands of low-mass protostars, 70% of which have disks (Hillenbrand et al. 1998). 
Some of the disks dbTQ clS close few tenths of a parsec from 6^ C Ori. 

Though evolving on substantially shorter timescales than low-mass stars, high mass 
stars still need a finite time to complete hydrogen and helium burning, and explode in a SN 
(Schaller et al. 1992; Romano et al. 2005). It takes ~ 40 Myr for a 8 M© B3 star and ~ 3 
Myr for an extremely massive 120 Mq 03 star to reach the SN stage respectively. A 40 M© 
star, such as 9^C Ori, evolves for ~ 5 Myr before it goes supernova (Figured]). 

The age of the ONC is estimated to be ^ 1 Myr (e.g. Hillenbrand 1997). The age of 
9^ C Ori is not precisely known, although this star is probably among the youngest members 
of the ONC, because early formation of 6^ C Ori would have halted star formation through 
rapid photoionization of the surrounding gas (O'Dell 2001; Boss & Goswami 2007). O'Dell 
(2001) argues that 6^C Ori is the youngest star in the ONC cluster. Palla & Stalher (2001) 
determined the age of another Trapezium star, 6^D Ori, and find an upper limit of 0.1 Myr. 
Here it is conservatively assumed that 6^ C Ori is 1 Myr old and will go SN 5 Myr after the 
onset of star formation in the ONC, i.e. ~ 4 Myr from now. 

The protoplanetary disks that are currently within a few tenths of a parsec from 6^C 
Ori will suffer severe mass-loss due to photoevaporation driven by the UV radiation from 
this and other massive stars in the region (Johnstone, Hollenbach & Bally 1998; Storzer 
& Hollenbach 1999; Richling Sz Yoke 2000; Balog et al. 2007). Typical mass loss rates 
due to photoevaporation are 10"'' M0/yr (Johnstone, Hollenbach & Bally 1998; Storzer & 
Hollenbach 1999). The typical lifetime of a minimum mass solar nebula of mass 0.01 Mq 
(Hayashi et al. 1985) in such an environment is therefore ~ 10^ yr. Johnstone, Hollenbach 
& Bally (1998) studied the specific case of 9^ C Ori and showed that, at a distance of 0.3 pc. 
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a protoplanetary disk would shrink to a size of about 1 AU in about 1 Myr, i.e. well before 
9"^ C Ori becomes a SN. Disk evaporation simply prevents injection of SN ejecta because disks 
in the vicinity of the SN have essentially disappeared before the SN explosion takes place. 

It is possible that disks adjacent to 6^ C Ori could partially survive photo-evaporation if 
their orbits around the massive star are highly elliptical rather than circular (e.g. Adams & 
Bloch 2005), or if the photo-evaporation rate decreases rapidly with time (e.g. Adams et al. 
2006). However, in 4 Myr from now, when 9^C Ori goes supernova, such surviving disks will 
be highly evolved, harboring large planetesimals as well as giant planets (e.g. Montmerle et 
al. 2006). In an Orion-like setting, planetary formation is accelerated by photoevaporation 
(Throop & Bally 2005), reducing accordingly the time window for SR injection. Disks which 
survive evaporation will therefore not receive SRs from the SN explosion until very late in 
their evolution, inconsistent with cosmochemical evidence for early delivery, which indicates 
that the injection happened within 1 Myr of disk formation (see section [2]). 

Over the next 4 million years, there will be few new low-mass stars forming within 1 pc of 
9^ C Ori. Star formation is certainly more vigorous during the first Myr of a molecular cloud 
lifetime and decreases sharply with time (see Lada & Lada 2003). After a few Myr, the gas in 
a molecular cloud has dissipated, and new star formation no longer takes place (Leisawitz, 
Bash & Thaddeus 1989; Elmegreen 2000; Hartmann, Ballesteros-Paredes & Bergin 2001; 
Lada & Lada 2003). Winds from massive stars and outflows from low-mass stars are primarily 
responsible for the quick dissipation of the molecular gas (Bally, Moeckel & Throop 2005). 
In fact, young clusters, such as NGC1333 and Serpens, contain relatively large populations 
of Class sources and are rich in bipolar flows, both features indicative of the earliest stages 
of star formation (e.g. Knee & Sandell 2000; Sandell & Knee 2001). In contrast, older 
clusters, such as IC348, contain few protostellar sources and outflows (Lada & Lada 2003). 
The fraction of disks in a given cluster, which is homologous to the star formation rate, also 
steeply decreases with the cluster age (see Figure 1 of Haisch, Lada & Lada 2001). 

The decrease with time of the star forming rate is even more dramatic in the immediate 
vicinity of a massive star (< few pc) than at the global molecular cloud scale. Hillenbrand 
(1997) argues that star formation in the Trapezium, i.e. in the vicinity of 9^C Ori, is now 
over. More generally, star formation cannot occur in the absence of molecular gas, which is 
ionized by the enhanced UV flux of massive stars. Wren & O'Dell (1995) estimate that the 
ionized HII region created by 9^C Ori is now 0.3 pc wide and is growing at the rate of 0.5 
km/s. Within the next 4 Myr, this region will expand to a ~ 2.3 pc wide HII cavity around 
9^C Ori, where star formation will be effectively halted. Hillenbrand & Hartmann (1998) 
estimate that the ~ 1 pc wide gaseous region surrounding 9^ C Ori will be photoevaporated 
on a timescale of ^ 1 Myr. 
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It is not only the ionizing power of massive stars, such as 6^ C Ori, that contributes to 
the sharp decrease in star formation rate in their immediate vicinity. The winds emanating 
from the massive stars create bubbles. The linear dimension of a wind bubble increases 
with the mass of the star (Chevalier 1999). A 20 Mq star creates a cavity as large as 11 pc 
in 7 Myr (Chevalier 1999), suggesting that the molecular gas surrounding 6^C Ori will be 
dissipated by winds in addition to photoionization on timescales significantly shorter than 
the lifetime of 6^C Ori. 

After a few Myr of evolution, star formation in ONC-like settings occurs mainly in 
photo-dissociation regions (PDRs) at the interface of the HII region and the molecular gas 
(Healy, Hester & Claussen 2004), a few parsecs away from the massive star. The fate of 
the ONC is well illustrated by the 2-3 Myr old NGC2244 cluster, whose most massive star, 
HD206267, has the same spectral type as 6^ C Ori. In NGC2244, star formation is occurring 
in the outskirts of the cluster, at distances 5-10 pc from the central 06 star (Hartmann 
2005). In addition. Reach et al. (2004) note that these stars, which formed ~ 2 Myr before 
HD206267 goes supernova, are the last generation of stars forming in the cluster. Low-mass 
stars formed 3 Myr after the onset of star formation in a given molecular cloud are too far 
away to be contaminated by SRs at any significant level. 

Figure [2] illustrates the discrepancy between the timescales of massive star evolution 
and low-mass star formation for an ONC-like setting. 



4. Probability of SR injection in a young disk by a nearby SN 

In this section, we calculate two numbers: (1) Within a cluster of a given size A^^,, what 
is the probability for a young disk (^ 1 Myr) attached to a low- mass star to be contaminated 
in SRs by a nearby supernova explosion at the level observed in the early solar system? (2) 
What is the probability for any protoplanetary disk in the Galaxy to be contaminated in SRs 
from a nearby SN explosion at the level observed in the solar system? The latter number 
will be calculated by integrating the former with the cluster size distribution (section 14.61) . 

The probability for a disk surrounding a low- mass star in a cluster of a given size (A^*) to 
be enriched in SRs by a nearby supernova at the observed solar system abundance without 
being disrupted by UV photoevaporation is given by the expression: 

-Psr(A^*) = / /YSo(tsN) /enrichment /survivors "^A/sn (^*)'^^SN, (3) 
J Ml 

where /yso is the fraction of young (^ 1 Myr) low mass stars (or Young Stellar Objects) 
present in the cluster at the time (tsN) of the supernova explosion, /enrichment is the fraction 
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of low-mass stars present in the enrichment zone, /survivors is the fraction of disks which 
have survived UV photoevaporation and VmsnI-^*) probabihty that the most massive 

star in a cluster containing A^* stars has a mass M^tss- M/ = 8 Mq and = 150 Mq are, 
respectively, the masses of the least and most massive star possibly going supernova (Kroupa 
& Weidner 2005). The dependence of /ysos, /enrichment, /survivors and Vmsn on and Msn 
will be described in the following sections. 



4.1. The most massive star in a cluster 

Using the Initial Mass Function (IMF), dN/dM oc M~^^^°'\ where N is the number of 
stars and M the star mass and with a = 1.5 (Scalo 1986), Williams & Gaidos (2007) show 
that the probability that the most massive star in a cluster containing A^* stars has a mass 
Msn is: ^ 

a/sNiV* f Mu\ " „-Ar(>MsN) 



where /sn = 0.003 is the total fraction of stars going supernova, i.e. having a mass larger 
than 8 Mq (Adams & Laughlin 2001). A/'(> Msn), the number of stars with mass larger 
than Msn, is obtained from the IMF and given by the expression: 

The dependence of VMsNi-^*) '^^ ^* shown in Figure [3] for two fiducial masses, Msn = 40 
Mq and Msn = 120 Mq. 



4.2. Fraction of stars in the enrichment zone 

Young stellar clusters are segregated in mass, with the more massive stars present at the 
center of the cluster (Hillenbrand & Hartmann 1998 and references therein). This core-halo 
structure appears to be primordial, based on morphological (Larson 1982) and dynamical 
(Hillenbrand & Hartmann 1998) evidence. It is thus reasonable to consider that clusters 
have a spherical symmetry, and that the massive star responsible for the putative injection 
of SRs is placed at the center of the cluster (see for example Rho et al. 2006). The number 
of stars in the enrichment zone is therefore given by the expression: 

/enrichment (A^*) = / p{r)ATcr'^dr / / p{r)AiTr'^dr, (6) 
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where p(r) is the stellar volume density and rmin and rmax define respectively the minimum 
and maximum acceptable distance between the disk and the SN, defining the enrichment 
zone (see section iV* and Rc(A^*) are the number of stars contained in the cluster and 
the cluster radius, respectively. 

Observations of ONC as well as IC348 show that the numbers of stars per unit area, 
i.e. the internal cluster density distribution, decreases approximately as 1/r (Hillenbrand & 
Hartmann 1998; Muench et al. 2003), in which case the stellar volume density, p{r), goes as 
1/r^. Inserting the latter expression into equation yields: 

f (M \ — ^max — ^min /„x 

Jenrichmentl,-' »*J o / Tvr \ ' V ' 

Using the embedded clusters compilation of Lada & Lada (2003), Adams et al. (2006) 
find that embedded clusters follow the law Rc ~ R300 a/-/V*/300, where Rc is the cluster size 
and R300 is an empirical parameter varying between 1 and 3 pc (Figure 2 of Adams et al. 
2006). This provides an estimate of the average surface density of clusters of different sizes, 
which can be formulated as: 

S = Njnil (8) 

We consider a wide range for E encompassing the range observed by Adams et al. (2006), 
i.e. S varies from 10 to 300 stars/pc^. Note that the ONC is characterized by Sqnc = 211 
stars/pc^, based on the cluster parameters (Rc = 2.06 pc and N,, = 2817 stars) of Hillenbrand 
& Hartmann (1998). Inserting equation ([8]) into equation ([7]) yields: 

J enrichment l^-' '* J / , , • ) 



4.3. Fraction of disks surviving photoevaporation 

Disks residing within a few tenths of a parsec of a massive O star will be UV photoe- 
vaporated on timescales significantly shorter than 1 Myr (e.g. Storzer Sz HoUenbach 1999). 
Because disks can have elliptical orbits around the massive star, and because disk-shrinking 
is limited by the decrease of the mass loss rate with the disk size (Adams et al. 2004), a 
non-zero fraction of disks might survive nearby massive stars. To estimate the fraction of 
disks that survive total photoevaporation we will turn again to Orion. 

Using the Hubble Space Telescope planetary camera, Johnstone, Hollenbach & Bally 
(1998) measured the size of 28 protoplanetary disks in the vicinity of 6^C Ori. Only nine of 
them have sizes larger than 35 AU, the estimated minimum size of the solar protoplanetary 
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disk (e.g. Gomes et al. 2005; Hartmann 2005). This means that, in the case of the ~ 1 
Myr old ONC, only a disk fraction /survivors = 0.32 has resisted to photoevaporation and are 
eligible to SN pollution. This number is compatible with the depletion of disks by a factor 
of ~ 3 observed by Balog et al. (2007) within 0.5 pc of O stars in the 2-3 Myr NGC 2244 
star forming region. 



4.4. Fraction of young low-mass stars at time t 



SN 

As discussed in section |2l cosmochemical constraints require the SN injection to occur 
early in the evolution of the protoplanetary disk. We conservatively estimated that injection 
had to occur in the first 1 Myr of disk evolution. The fraction of low-mass stars younger 
than = 1 Myr at time t is: 

/Yso(t) = /o Ht')dt' I ^it')dt' for t ^ e 
/Yso(t) = lieHm' / J.^i'indt' fort^e (10) 

where iplt) is the star forming rate. 

Star forming rates are not very well constrained and are the subject of debate among 
astronomers. We therefore consider three different expressions for ?/'(t), a combination that 
is likely to bracket reality. 

The first expression for ^(t) is a step function with a total star formation timescale T 
= 4 Myr: 

' ^(t) = ipi for t ^ T 
^{t) = for t ^ T ^^^'> 

Expression f|TT]) gives an overestimated measure of the star formation rate, because star 
forming regions are known to be more active during their first million years than during 
subsequent times (e.g. Hartmann, Ballesteros-Paredes & Bergin 2001). 

The second expression for iplt) is a linearly decreasing function, terminating on a 
timescale of T = 4 Myr, due to molecular gas dissipation (Leisawitz, Bash & Thaddeus 
1989; Elmegreen 2000; Hartmann, Ballesteros-Paredes & Bergin 2001; Lada & Lada 2003). 

V'(t) = ^2 X (1 - I) for t ^ T 
ip{t) = fort^T 



The third expression for ?/'(t) is an exponentially decreasing function with a character- 
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istic star formation timescale of r = 0.56 Myr: 

tP{t) = ^0 X (13) 

Expression (JT3ll is determined from a fit of star formation rates vs. ages for a suite of observed 
clusters shown in Figure HI We emphasize that all the star forming regions considered in 
this plot have similar sizes (i.e. similar contents of gas). This expression describes a rather 
extreme, though possible, star formation history with most stars being formed during the 
first million year of evolution. 

Our preferred model corresponds to the linearly decreasing law for z/'(t) as it describes 
an intermediary situation between the two extreme cases, the step function and the expo- 
nentially decreasing law. 

The fraction of low- mass stars (calculated with equation fllOj) ) present in the cluster as 
a function of time, /Yso(t), is shown in Figure [5] for the three different star formation rates 
considered. The sharp cut-offs at ~ 5 Myr in Figure [5] for the linear and the step function 
star formation rates are due to finite duration of star formation in an embedded cluster 
(Lada & Lada 2003 and references therein). For these cases, there are no low-mass stars 
young enough to be enriched in the cluster at a time t = T -|- 6 = 5 Myr. 

For disk contamination to occur, the relevant quantity is the fraction of low-mass stars 
at the time tsN of the supernova explosion, i.e. /Yso(tsN)- Evolutionary models for massive 
stars (Schaller et al. 1992; Romano et al. 2005) provide the time tsN at which a massive 
star goes supernova through the expression: 

log(tsN) = 1.4/(logMsN)'■^ (14) 
where tsN and Msn are expressed in Myr and in solar masses respectively. 



4.5. Injection probability in a cluster containing A^^, stars: Results 

Putting together the expressions defined above, we obtain: 

Psr{N.) = /survivors /YSo(tsN) ^^'"'^^ ~ ^'^^ ^ VmsA^^MsN, (15) 

J Ml V -'V* 

where /Yso(tsN) and "PmsnI^*) are given by equations (|TOl) . IHM and (jlj) respectively. With 
S = SoNC = 211 stars/pc^ (section and r^i^ = 0.1 pc and r^ax = 0.4 pc (section [2]), 
we obtain a set of curves for Psr(A^^,) shown in Figure [61 Our curves are slightly below that 
calculated by Williams & Gaidos (2007) because we take explicitly into account the disk 
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photoevaporation, and because we model the timing of injection in a fashion dictated by the 
cosmochemical constraints, i.e. we consider that injection happened during the first Myr of 
disk evolution rather than at any time (section 2). 

The probability for SR enrichment peaks at ~ 10 000 stars. It is however very low as 
it reaches a maximum value of ~ 0.4 % for the most favorable case, corresponding to a step 
function law for the star formation rate ip{t). It peaks at ~ 0.03 % for the exponential form 
of the star formation rate, which is based on the data shown in Figure HI Our preferred model 
(a linearly decreasing function for ip{t)) gives an intermediary result with a maximum value 
of ~ 0.3 %. Increasing E to its maximum value, E = 300 stars/pc^, does not substantially 
change these results, because Psr(-^*) varies smoothly with the square root of E (see equation 
ffT^ ) and all our calculations are normalized to the high stellar density of the ONC. 



To calculate the probability of SR injection for any protoplanetary disk attached to a 
low-mass star in the Galaxy, Vqal, we will integrate Psr(^*) over the probability for a star 
to form in a cluster of size iV*, which scales as 1/A^* (Elmegreen & Efremov 1997; Adams & 
Laughlin 2001). In other words. 



where Nmin = 100 and N^ax = 5 x 10^ are the minimal and maximal sizes of stellar clus- 
ters respectively (McKee & Williams 1997; Williams & Gaidos 2007). The normalization 
constant, /C = 9.39 x 10~^, is calculated using the observation that ~ 80 % of stars form in 
clusters larger than 100 members (Lada & Lada 2003). Figure [7] shows Vqal for the three 
different star forming rates discussed in section 14.41 In our preferred case, where E = Eqnc 
= 211 stars/pc^ and ijj{t) follows a linearly decreasing law, the probability for any disk to 
receive SRs from a nearby SN at the level observed in the solar system is 1 x 10~^ (grey 
circle in Figure [7]) . This probability increases with the square root of E as expected from 
equation fllSp . In the most favorable case, i.e. adopting a step function for the star formation 
rate and with E = 300 stars/pc^, the probability for any diskS to receive SRs from a nearby 
SN at the level observed in the solar system is 3 x 10"'^. 



4.6. Probability of injection for any disk 




(16) 
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4.7. Multiple supernovae 



If the cluster size is large enough, multiple supernovae can occur within the lifetime of 
the cluster. Stars with mass below 40 Mq can not contribute to the inventory of SRs because 
they need too long timescales to evolve (see section [3]). Multiple supernovae injection of SRs 
therefore occur only in clusters large enough to contain at least one star whose mass is larger 
than 40 M0. This condition is satisfied only for clusters containing more than iV* = 10^ 
stars (Equation ([5]) and Figure [8]). 

To calculate the maximum contribution of nearby massive stars to a young protoplane- 
tary disk, we assume that, for clusters containing more than = 10^ stars, every low-mass 
star will witness a supernova explosion in its youth (^ 1 Myr), i.e. that: 



Introducing this condition in equation ( ITSl) yields a stringent upper limit, as this statement 
is probably true only for the most massive clusters (N,, ^ 10^), where there are enough 
stars to include tens of supernovae whose progenitors stars have masses larger than 40 Mq. 
With this generous permission, we calculate the probability for any young low-mass star 
protoplanetary disk to be contaminated in SRs by a nearby supernova explosion to be 6 x 
10~^ in the most favorable case (^(t) described by a step function and E = 300 stars/pc^). 



The probability of a protoplanetary disk to receive SRs from a nearby SN explosion 
would be higher if one considered a star forming region more massive than the ONC. Such a 
region would contain more massive stars than the ONC according to the initial mass function 
(Salpeter 1955). For example, a 120 Mq star would be ready to go SN and deliver SRs to 
a nearby disk after only ~ 3 Myr of evolution (Figure [1] and equation (fT4|) ). But, 3 Myr 
after the onset of star formation, the star formation rate has already decreased significantly 
as can be appreciated in Figure |H Additionally, extremely massive molecular clouds are 
rare, at least in the vicinity of the Sun. Within 2 kpc of the Sun, where the number of 
embedded clusters is large enough to make statistically robust observations, the ONC is the 
most massive example (Lada & Lada 2003). Finally, a more massive cluster contains more 
massive stars which are more efficient in dissipating molecular gas by ionization and photo- 
dissociation (e.g. Bally, Moeckel & Throop 2005). Despite these difficulties, the massive 




1. 
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5. 



Discussion 



5.1. Beyond Orion: The Carina Nebula 
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Carina Nebula star-forming region is now proposed by some authors (e.g. Smith & Brooks 
2007), instead of Orion, as a good analog for the astrophysical environment of our solar 
system formation. 

The Carina Nebula (NGC3372) is a star forming region situated at 2.3 kpc from the 
Sun and containing ~ 65 O stars (Smith 2006). rj Carinae, the most massive star contained 
in the Carina Nebula (^ 100 M©, Davidson & Humphreys 1997), needs ~ 3 Myr to go 
supernova (Smith 2006). Since the age of the cluster is ~ 3 Myr, assuming it formed at the 
beginning of the cluster's life, rj Carinae can go supernova anytime now (Smith 2006). Most 
of the molecular gas in the vicinity of rj Carinae has however been cleared away, creating an 
extended HII region where star formation is halted (Smith et al. 2000 and Figure 5 of de 
Graauw et al. 1981). 

A few evaporating protoplanetary disks were tentatively detected in the Carina Nebula 
and await confirmation (Smith, Bally & Morse 2003). From the coordinates given by Smith, 
Bally & Morse (2003), we calculate that the protoplanetary disk closest to rj Carinae lies at 
a minimum distance of 2.4 pc, a factor of ~ 10 too large compared to the distance needed to 
inject SRs in sufficient number. In addition, the putative protoplanetary disks detected by 
Smith, Bally & Morse (2003) are significantly more massive (by a factor of ~ 100) than the 
protoplanetary disks observed in Orion (Johnstone, HoUenbach & Bally 1998), or than the 
expected solar protoplanetary disk (0.013 Mq, Hayashi et al. 1985), suggesting that they 
belong to young intermediate-mass stars rather than to solar-like stars. 

The scarcity of solar-like protoplanetary disks in the vicinity of rj Carinae is due to the 
extreme ionizing power of that star which emits Q{H) ~ 10^°-^ ionizing photons per second, 
which can be compared to the Q{H) ~ 10^^ photons per second emitted by 6^C Ori (O'Dell 
2001). It confirms that in the immediate vicinity of very massive stars, protoplanetary disks 
are extremely rare because (1) those accompanying low-mass stars which formed contem- 
poraneously with massive stars were photoevaporated or formed planets (Throop & Bally 
2005) and (2) the efficient gas clearing by winds and UV emission has halted new low-mass 
star formation. 

Triggered star formation in the Carina Nebula occurs in the gas-rich region named 
South Pillars (Smith & Brooks 2007). Smith & Brooks (2007) argue that young stars in 
the South Pillars region will be pelted by supernova ejecta containing SRs. While it is true 
that these young stars will be exposed to the supernova ejecta, the South Pillars lie 20 pc 
away from rj Carinae (Smith et al. 2000), almost two orders of magnitude further away 
than the distance required to incorporate SRs at the abundance observed in the solar system 
(see section |2]). Massive stars which form now in the South Pillars will need several Myr 
to explode and deliver SRs to the interstellar medium and therefore suffer from the same 
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difficulties discussed above. 



5.2. On our probability estimate 

In calculating the different probability estimates for the injection of SRs into a young 
protoplanetary disk attached to a low-mass star, we made several assumptions and simplifi- 
cations that we discuss in the following. 

5.2.1. The enrichment zone 

We used a minimum distance between the supernova and the disk, rmin = 0.1 pc, pro- 
posed by Ouellette, Desch & Hester (2007) from models simulating the interaction of a 25 
Mq SN with a protoplanetary disk. A more massive star, more likely to achieve that goal, 
would have a more massive ejecta and therefore a more disruptive effect (through momen- 
tum transfer) on the disk (Chevalier 2000). The minimum SN-disk distance we adopted is 
therefore a lower limit, which turns into an upper hmit for /enrichment (see equation [3]) and 
therefore into an upper limit for the probability of SR delivery by a nearby SN in a young 
disk. 

The distance r varies with the square root of the injection efficiency {rf) and with the 
inverse square root of the initial abundance of short-lived radionuclides (Mgs) (equation 
([2])). The initial abundance of SRs is known within a factor of a few (Gounelle 2006). 
Therefore, changes in their adopted abundance will not change much the results. In the 
case of injection into a protoplanetary disk, Ouellette, Desch & Hester (2007) show that the 
injection efficiency is below 1 % for the gas fraction of the ejecta. Though they suggest that 
SRs could be injected in the disk in the form of dust grains, it remains to be demonstrated 
how plausible dust injection is, especially given that SN dust condensation efficiency (defined 
as the ratio of mass of refractory elements condensed into dust to that of refractory elements 
in the ejecta) is ^ 0.12 (Sugerman et al. 2006 and references therein). Ercolano et al. (2007) 
estimate an upper limit as low as 4 x 10~^ for the dust condensation efficiency of supernova 
SN 1987A. If Tj is reduced to one order of magnitude in equation ([2]), it will reduce the 
maximum acceptable distance by a factor of vTo ~ 3, i.e. from 0.4 pc to ~ 0.15 pc. The 
maximum distance rmax calculated in section [2] is thereore a strict upper limit, which turns 
into an upper limit for /enrichment (see equation [3]) and, accordingly, for the probability of SR 
delivery by a nearby SN in a young disk. 



-17- 



5.2.2. Disk UV photoevaporation 

To estimate the fraction of disks which survived UV photoevaporation(/survivors), we 
calculated from the data of Johnstone, Bally & Hollenbach (1998) the fraction of disks in 
Orion which are larger than the minimum inferred original size of the solar protoplanetary 
disk (35 AU, Gomes et al. 2005). This provides a strong upper limit for /survivors for three 
reasons. First, the ONC is only ~ 1 Myr old and at the time of 6^ C Ori explosion (4 Myr 
from now) many more disks will have their sizes reduced below 35 AU. Second, this estimate 
does not include the disks which have already fully evaporated and will not capture any SRs. 
Third, more massive stars, more likely to deliver SRs, have also larger UV fluxes which will 
result in a more rapid photoevaporation (Johnstone, Bally & Hollenbach 1998). The fraction 
of disks that survive UV photoevaporation, /survivors, which we used in our calculations is 
therefore an upper limit, as is accordingly the probability of SR delivery by a nearby SN in 
a young disk. 

5.2.3. Temporal coincidence of disks and supernovae in a molecular cloud 

In calculating /yso; we considered ipit), the temporal evolution of the star formation 
rate, at the global molecular cloud scale, and did not take into account the dissipation of 
molecular gas in the immediate vicinity of the most massive star. This is a simplification 
as discussed in section [31 Because massive stars create very rapidly after their formation a 
HII region in which star formation is impossible, star formation in molecular clouds occurs 
in the outskirts of HII regions, at several parsecs from the massive star (Reach et al. 2004; 
Balog et al. 2007; Hartmann 2005). In other words, after 2-3 Myr of evolution, the star 
formation rate, ip{t), in the immediate vicinity (< 2 pc) of any massive star is virtually zero. 
By assuming non-zero star formation rates, we definitely calculated a strong upper limit for 
/yso and therefore for the probability of contamination of a young protoplanetary disk by a 
nearby supernova. 

In calculating Psr(A^*) (equation ([3])) and /yso (equations (|T0|) and f|T^ ). we implicitly 
assumed that star formation was coeval in the considered star forming region. In fact, it 
is only if the most massive stars in a cluster formed in advance relative to low-mass stars 
that the probability of injection can increase. In such a case, the rapid evolution timescale 
of protoplanetary disks could be reconciled with the slower evolution timescale of massive 
stars (Figure [2j) and the destructive effect of massive stars would play a less important role. 

It is commonly argued that molecular clouds have short lives (~ 4 Myr) and that star 
formation within them proceed as soon as the cloud is assembled (e.g. Ballesteros-Paredes 
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et al. 2007). This implies that the star formation age spread is small (though non zero, e.g. 
Hartmann 2001). On the other hand Palla & Stalher (2000) have argued that molecular 
clouds live some tens of Myr and that some stars form in advance compared to the majority 
of stars. Even if this is correct, it would not solve the timescale problem, because all stars 
which possibly formed early relative to the majority of stars in the cluster have (at least in 
the case of the ONC) masses below 0.3 Mq (Palla et al. 2007). Such stars are far from being 
massive enough to go SN and deliver SRs to a nearby protoplanetary disk. We note that, 
in general, high-mass stars are likely to form last in a cluster (Bally et al. 1998; Storzer & 
Hollenbach 1999; O'Dell 2001; Kroupa & Weidner 2005; Henriksen 1986; Boss & Goswami 
2007). This implies again that the probability estimates given in the previous sections are 
upper limits. 



5.2.4- Multiple supernovae 

The contamination probability in the case of multiple supernovae was explicitly overes- 
timated in assuming (equation f|T71) ) that every low-mass star in the cluster will witness a 
supernova explosion in its youth (T ^ 1 Myr). The overall contamination probability was 
also implicitly overestimated by keeping /enrichment identical to that calculated in the single 
supernova case. This is because, in calculating /enrichment (section we assumed that the 
most massive star is located at the center of the cluster, where the density of stars is highest. 
In the case of multiple supernovae, the additional massive stars are located in the cluster 
periphery (as only the most massive star occupies the cluster center) where low-mass stars 
are rare, decreasing accordingly /enrichment- 

Even with those gross approximations, the effect of multiple supernovae is to increase 
the probability by a factor of only a few. This factor of a few would cancel if the physical 
effect of many more massive stars were properly taken into account. More massive stars 
would lead to increased photoevaporation power, clearing very large HII regions where star 
formation is halted. For example the HII region surrounding rj Carinae has an expansion 
radius of 110 pc (Smith et al. 2000). In addition, the presence of many more massive 
stars would also increase the probability of close stellar encounters, and therefore planetary 
disruption (see below). For all these reasons, we will ignore the increase by a factor of a few 
in the calculated probability due to multiple supernovae. 
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5.2.5. Planetary disruption 

A low-mass star born in a large cluster where many massive stars are present can endure 
close stellar encounters, which have the possibility to disrupt the orbits of planets and small 
bodies. This was modeled by Adams & Laughlin (2001). They showed that for clusters 
larger than ~ 2500 members, the cluster density is such that stellar encounters would modify 
the orbits of the giant planets and the Kuiper Belt. This led Adams & Laughlin (2001) to 
suggest that the Sun formed in a cluster smaller than a few thousand stars. In doing so, they 
considered relatively long relaxation timescales for clusters. This was criticized by Hester 
& Desch (2005). Though small clusters live only up to 10 Myr before disruption, there is a 
non zero fraction of clusters, corresponding to the most massive ones, which live hundreds 
of Myr (Lada & Lada 2003), similarly to the lifetimes of the Pleiades. Kroupa et al. (2001) 
suggested on the basis of a dynamical study that the ONC might become an open cluster 
similar to the Pleiades, reinforcing the idea that for large clusters the relaxation timescales 
adopted by Adams & Laughlin (2001) is correct. Because other complications might arise 
(e.g. Megeath et al. 2007) we did not, however, take explicitly into account the effect of 
close stellar encounters in our calculations. Neglecting these effects for large clusters implies 
that the probability of SR delivery by a nearby SN to a young disk, which we calculated 
above, is an upper limit. 



6. Implications for the origin of short-lived radionuclides 

6.1. The origin of short-lived radionuclides 

The probability for any protoplanetary disk attached to a young low-mass star to receive 
SRs from a nearby SN at the level observed in the solar system is ~ 1 x 10"'^ in our preferred 
case. As shown in section 15.21 this value is a stringent upper limit. Although improbable 
does not mean impossible, this very low probability suggests that other sources should be 
considered for explaining the overabundance (compared to the interstellar medium average 
value) of SRs in the early solar system (Wadhwa et al. 2007). Some of the SRs, such as 
^^Al, ^^Cl, ^^Ca and ^^Mn can be made by in situ irradiation together with ^Be and ^°Be 
(e.g. Gounelle et al. 2006), and therefore do not necessitate further elucidation. Iron-60, 
on the other hand cannot be made by local irradiation (Lee et al. 1998). It is therefore 
necessary to find a source of ®°Fe, within a plausible astrophysical context. Besides the 
nearby supernova scenario, which seems to be unlikely, there are two remaining possibilities 
for the ^'^Fe presence in the solar system: (1) a distant SN or (2) an inherited origin resulting 
from the contributions of many SN. 
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A distant supernova (i.e. at a few parsecs from the solar system) was invoked in the 
past as the source of ^°Fe and other short-hved radionuchdes (Cameron & Truran 1977; 
Cameron et aL 1995). In that context, it was also assumed that the supernova shockwave 
triggered the collapse of the dense core. This proposition has a number of problems. First, 
very special conditions are needed to inject short-lived radionuclides in a dense molecular 
cloud core without disrupting it (e.g. Boss & Vanhala 2000 and references therein). The 
supernova shockwave needs to impact the core with a fine-tuned velocity of ~ 20 km/s. 
A higher shockwave velocity would disrupt the dense core, while a slower one would fail 
to induce collapse and inject radioactivities. Second, injection calculations do not take into 
account the filamentary nature of molecular clouds, which consist of the fractal juxtaposition 
of high density clumps and low density interclumps matter (e.g. Rho et al. 2006), but treat 
the cloud as an homogeneous and dense matter. A supernova shockwave within a molecular 
cloud will follow the path of less resistance, and will probably avoid the denser regions and 
instead expand into the interclump matter (Chevalier 1999) delivering very few, if any at 
all, short-lived radionuclides to molecular cloud cores. Third, though the larger size of a 
molecular cloud core compared to a protoplanetary disk allows it to collect SRs from a SN 
ejecta at larger distances (see equation [2]), it cannot be much further away than a few pc 
(Looney, Tobin & Fields 2006), and many of the difficulties described above hold in that 
case too. 



6.2. The solar system initial content of ^°Fe 

For some years, a high inferred initial ^^Fe/^^Fe ratio (e.g. Mostefaoui et al. 2005) was 
considered as a smoking gun for ruling out an inherited origin for ^°Fe (Hester & Desch 
2005). The experimental situation has now changed as was recently discussed at the Hawaii 
meeting on the chronology of meteorite^. The maximum solar system ^°Fe/^^Fe initial ratio 
of 1 X 10~^ (e.g. Wadwha et al. 2007) is probably too high an estimate. Below we give a 
summary of the various estimates for the initial ^°Fe/^^Fe ratio in the solar system. 

Measurements of Ni isotopes in CAIs are challenging because CAIs are rich in Ni nucle- 
osynthetic anomalies, which can blur the effect of ^°Fe decay (Birck 2004). An upper limit 
on the ^°Fe/^^Fe ratio of 1.6 x 10^^ was given by Birck & Lugmair (1988) for an AUende 
CAI with a Fe/Ni ratio of ~ 20. We note that no isochron was reported for this object, 
and that the authors expressed caution about the ^°Fe decay origin of the small (~ 0.01 %) 
anomaly detected. A two point internal isochron for an AUende CAI measured by Quitte et 



http://www.lpi.usra.edu/nieetings/metchron2007 
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al. (2007) gave a lower limit for the initial ^^Fe/^^Fe ratio of 3 x 10^^ 

Measurements on the silicate portion of chondrules from primitive meteorites, including 
Semarkona, yield ^^Fe/^^Fe ratios between 1.7 and 3.2 x 10"'^ (Tachibana et al. 2006; 
Tachibana et al. 2007), while a measurement from a Semarkona troilite (FeS) gives ™Fe/^^Fe 
~ 9 X 10-'^ (Mostefaoui et al. 2005). High ^^Fe/^^Fe initial ratios in troilite can be due to 

Fe-Ni redistribution in the sulfides during later alteration processes (Guan et al. 2007). It is 
indeed difficult to understand how troilite would have formed with a higher initial ^°Fe/^^Fe 
ratio than silicates in chondrules from the same meteorite. If an hypothetic time delay of 
1.6 Myr is assumed between the formation of the solar system and the closure of the Fe-Ni 
system in ordinary chondrites (Connelly et al. 2007), initial ^^Fe/^^Fe ratios between 3.6 
and 6.7 x 10~^ can be calculated from the measurements of Tachibana et al. (2006) and 
Tachibana et al. (2007). 

Differentiated meteorites add to the complexity of the Fe-Ni system. The absence of 
^°Fe evidence in minerals from eucrites and angrites with high Fe/Ni ratios lead Sugiura et 
al. (2006) to put an upper hmit of 1.5 x 10"''' on the initial ^°Fe/^^Fe ratio. Quitte et al. 
(2006) reported values of (1.3 ± 0.8) x 10"^ and (8.2 ± 2.6) x lO"'^ for the ^^Fe/^^Fe at 
the time of formation of the angrites Sahara-99555 and d'Orbigny. Combined with Pb-Pb 
ages of ~ 3.3 and ~ 2.6 Myr after CAIs respectively (Amelin et al. 2006; Amelin 2007), 
these estimates lead to an upper value for the initial ^"^Fe/^^Fe of 1 x 10^^. The study of 
the eucrites Chervony Kut and Juvinas lead Shukolyukov & Lugmair (1993) to put an upper 
hmit of 2 x lO""^ for the initial ^^Fe/^^Fe ratio. 

Although there is at present no consensus on the initial value of ^°Fe in the early solar 
system, a low value for the ^°Fe/^^Fe ratio of a few times 10~^ seems preferred. 



6.3. An inherited origin for ^°Fe 

While a high (1 x 10~^) initial ^°Fe/^^Fe ratio could be interpreted as a strong evidence 
for a nearby SN contamination, revision to a lower initial value (^°Fe/^^Fe ~ 3 x 10~^) 
changes the outcome because it is compatible with a galactic background origin for ^°Fe. 

The present average galactic value of the ^^Al/^^Al ratio has been estimated to be 8.4 
X 10~^ (Diehl et al. 2006). Recent measurements of 7-ray lines in the interstellar medium 
with the Integral satellite have yielded a ratio ''°Fe/^''Al = 0.148 (Wang et al. 2007). Using 
an elemental ratio ^^Al/^^Fe = 0.109 (Lodders 2003) gives an average galactic value for the 
60pg^56pg of 1.4 X 10~^. This latter value is identical to the lower estimates of the early 
solar system inferred from meteorites (see above). Even if the initial solar system ^^Fe/^^Fe 
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ratio was as high as 4 x 10 ^ (e.g. Quitte et al. 2007), fluctuations around the average value 
could easily explain the presence of ^°Fe due to an inherited origin in the solar system. 

It is important to note that this estimate is that of the present interstellar medium. It 
might have been significantly different 4.6 Gyr ago. In addition, this number is averaged on 
the entire inner Galaxy. Heterogeneities in the ^^Al abundance in the Galaxy (Knodlseder 
et al. 1999) suggest that both ^^Al and ^°Fe abundances can fluctuate around the average 
value calculated above. Such calculations show that if the initial solar system ^°Fe/^^Fe ratio 
was only a few times 10~^, the need for a single supernova vanishes and an inherited origin 
can account for the ^°Fe found in the early solar system. Though ^^Al decays faster than 
^°Fe, an inherited origin cannot be excluded for this radionuclide either!^. 

We are especially grateful to Fred Adams and Thierry Montmerle for their generous and 
stimulating inputs throughout the writing of that paper. We thank Edouard Audit, Jerome 
Bouvier, Roland Diehl, Jean Duprat, Eric Feigelson, Lee Hartmann, Patrick Hennebelle, 
Shu-Ichiro Inutsuka, Doug Johnstone, Mordecai Mac Low, Kevin McKeegan, Alessandro 
Morbidelli, Manuel Petitat, Ghylaine Quitte, Frangois Robert, Sienny Shang, Frank Shu, 
Vincent Tatischeff and Edward Young for illuminating discussions and sharp criticisms. An 
anonymous reviewer helped to improve the manuscript. This study made use of the ADS 
abstract service and was partially funded by the Programme National de Planetologie (PNP) 
and the CNRS fund France-Etats-Unis. 

REFERENCES 

Adams, F. C. & Bloch, A.M. 2005, ApJ, 629, 204 
Adams, F. C. & Laughlin, G. 2001, Icarus, 150, 151 

Adams, F. C., Hollenbach, D., Laughlin, G. & Gorti, U. 2004, ApJ, 611, 360 
Adams, F. C., Proszkow, E. M., Fatuzzo, M., & Myers, P. C. 2006, ApJ, 641, 504 
Amelin, Y., Wadhwa, M. & Lugmair, G. W. 2006, LPSC, 37, #1790 

^The recent claim for a late injection of ^"Fe (Bizzarro et al. 2007) was based on ^°Ni deficits (relative 
to chondrites) in iron meteorites and pallasites, which failed to be confirmed by other groups (Andrews et 
al. 2007; Regelous et al. 2007; Dauphas et al. 2008). There is at present no evidence for late injection of 
iron-60. 



-23- 



Amelin, Y. 2007, LPSC, 38, #1338 

Andrews, K. T., Schiller, M., Bizzarro, M., & Baker, J. A. 2007, Workshop on the Chronol- 
ogy of Meteorites and the Early Solar System, held November 5-7, 2007 in Kauai, 
Hawaii. LPI Contribution No. 1374, p.24-25, 24 

Ballesteros-Paredes, J. et al. 2007 in Protostars and Planets V, ed B.Reipurth, D. Jewitt & 
K. Keil (Tucson: University of Arizona Press), 63 

Bally, J., Moeckcl, N. & Throop, H. 2005, in Chondrites and the Protoplanetary Disk, ed A. 
N. Krot, E. R. D. Scott & B. Reipurth (San Francisco: ASP Conference Series), 81 

Bally, J., Testi, L., Sargent, A., & Carlstrom, J. 1998, AJ, 116, 854 

Balog, Z. et al. 2007, ApJ, 660, 1532 

Birck, J.-L. 2004, in Geochemistry of non-traditional stable isotopes, ed C. M. Johnson, B. 
L. Beard & F. Albarede (Washington D.C.: Mineralogical Society of America), 25 

Birck, J. L. & Lugmair, G. W. 1988, EPSL, 90, 131 

Bizzarro, M., Ulfbeck, D., Trinquier, A., Thrane, K., Connelly, J. N., & Meyer, B. S. 2007, 
Science, 316, 1178 

Boss, A. P. & Goswami, J. 2007, in Meteorites and early solar system II, ed D. S. Lauretta 
& H. Y. McSween (Tucson: University of Arizona Press), 171 

Boss, A. P. & Vanhala, H. A. T. 2000, SSR. 92, 13 

Bottke, W. F., Nesvorny, D., Grimm, R. E., MorbideUi, A. & O'Brien, D. P. 2006, Nature, 
439, 821 

Bouvier, A., Bhchert-Toft, J., Moynier, F., Vervoort, J. D., & Albarede, F. 2007, 
Geochim. Cosmochim. Acta, 71, 1583 

Busso, M., Gallino, R. & Wasserburg, G. J. 2003, Publications of the Astronomical Society 
of Australia, 20, 356 

Cameron, A. G. W. & Truran, J. W. 1977, Icarus, 30, 447 

Cameron, A. G. W., Hofiich, P., Myers, P. C. & Clayton, D. D. 1995, ApJ, 447, L53 
Chevalier, R.A. 1999, ApJ, 511, 798 
Chevalier, R. A. 2000, AJ, 538, L151 



-24- 



Connelly, J. N., Amelin, Y., Krot, A. N., & Bizzarro, M. 2007, Workshop on the Chronol- 
ogy of Meteorites and the Early Solar System, held November 5-7, 2007 in Kauai, 
Hawaii. LPI Contribution No. 1374, p.46-47, 46 

Dauphas, N., et al. 2008, ApJ, Accepted, 

Davidson, K., & Humphreys, R. M. 1997, ARA&A, 35, 1 

de Graauw T. et al. 1981, AA, 102, 257 

Diehl, R., et al. 2006, Nature, 439, 45 

Elmegreen, B. C, & Efremov, Y. N. 1997, ApJ, 480, 235 

Elmegreen, E.G. 2000, ApJ, 277 

Ercolano, B., Barlow, M. J. & Sugerman, B. E. K. 2007, MNRAS, 375, 753 
Galy, A., Young, E. D., Ash, R. D. & O'Nions, R. K. 2000, Science, 290, 1751 
Ghosh, A. & McSween Jr, H. Y. 1998, Icarus, 134, 187 

Gomes, R., Levison, H. P., Tsiganis, K. & MorbideUi, A. 2005, Nature, 435, 466 

Goswami, J. N., Marhas, K. K., Chaussidon, M., Gounelle, M., & Meyer, B. S. 2005, in 
Chondrites and the Protoplanetary Disk, ed A. N. Krot, E. R. D. Scott & B. Reipurth 
(San Prancisco: ASP Conference Series), 485 

Gounelle, M. & Meibom, A. 2007, ApJ, 664, L123 

Gounelle, M. 2006, NAR, 50, 596 

Gounelle, M., Shu, P. H., Shang, H., Glassgold, A. E., Rehm, K. E. & Lee, T. 2006, ApJ, 
640, 1163 

Gounelle, M. & Russell, S. S. 2005a, CCA, 69, 3129 

Gounelle, M. & Russell, S. S. 2005b, in Chondrites and the protoplanetary disk, ed A. N. 
Krot, E. R. D. Scott & B. Reipurth (San Prancisco: ASP Conference Series), 548 

Guan, Y., Huss, G. R. & Leshin, L. A. 2007, CCA, 71, 4082 

Haisch, K.E., Lada, E. A., Lada C.J. 2001, ApJ, 553, L153 

Hartmann, L.W. 2001, AJ, 121, 1030 



-25- 

Hartmann, L.W., Ballesteros-Paredes, J. & Bergin, E. A. 2001, ApJ, 562, 852 

Hartmann, L.W. 2005, in Chondrites and the Protoplanetary Disk, ed A. N. Krot, E. R. D. 
Scott & B. Reipurth (San Francisco: ASP Conference Series), 131 

Hayashi, C, Nakazawa, K. & Nakagawa, Y. 1985, in Protostars and planets II, ed D. C. 
Black & M. S. Matthews (Tucson: University of Arizona Press), 1100 

Healy, K.R., Hester, J.J., Claussen, M.J. 2004, ApJ, 610, 835 

Henriksen, R.M. 1986 , ApJ, 310, 189 

Hester, J. J. & Desch, S. J. 2005, in Chondrites and the protoplanetary disk, ed A. N. Krot, 
E. R. D. Scott & B. Reipurth (San Francisco: ASP Conference Series), 107 

Hillenbrand, L. 1997, A J, 113, 1133 

Hillenbrand, L. & Hartmann, L.W. 1998, ApJ, 492, 540 

Hillenbrand, L. et al. 1998, AJ, 116, 1816 

Kastner, J. H. & Myers, P. C. 1994, ApJ, 421, 605 

Knee, L. B. C, & Sandell, G. 2000, A&A, 361, 671 

Kn5dlseder, J., et al. 1999, A&A, 345, 813 

Kroupa, P., & Weidner, C. 2005, Massive Star Birth: A Crossroads of Astrophysics, 227, 
423 

Kroupa, P., Aarseth, S., & Hurley, J. 2001, MNRAS, 321, 699 
Johnstone, D., HoUenbach, D. & Bally, J. 1998, ApJ, 499, 758 
Lada, C. J. & Lada, E. A. 2003, ARAA, 41, 57 
Larson, R. B. 1982, MNRAS, 200, 159 

Lee, T., Shu, F. H., Shang, H., Glassgold, A. E. & Rehm, K.E. 1998, ApJ, 506, 898 
Leisawitz, D., Bash, F.N., Thaddeus, P. 1989, ApJS, 70, 731 
Lodders, K. 2003, ApJ, 591, 1220 

Looney, L. W., Tobin, J. J., & Fields, B. D. 2006, ApJ, 652, 1755 



- 26 - 

Luhman, K.L. 2004, ApJ, 602, 816 

Markowski, A., Leya, I., Quitte, G., Amnion, K., Halliday, A. N. & Wieler, R. 2006, EPSL, 
250, 104 

McKee, C. F., & Williams, J. P. 1997, ApJ, 476, 144 

McKeegan, K. D. & Davis, A. M. 2004, in Treatise on Geochemistry, ed H. D. Holland & K. 
K. Turekian (Oxford: Elsevier- Pergamon), 431 

Megeath, S. T., Gaidos, E., Hester, J. J., Adams, F. C., Bally, J., Lee, J. -., & Wolk, S. 
2007, ArXiv e-prints, 704, arXiv:0704.1045 

Meyer, B. S. & Clayton, D. D. 2000, SSR, 92, 133 

Montmerle, T., Augereau, J. C., Chaussidon, M., Gounelle, M., Marty, B. & Morbidelli, A. 
2006, EMP, 98, 39 

Mostefaoui, S., Lugmair, G. W. & Hoppe, P. 2005, ApJ, 625, 271 

Muench, A.A., Lada, C.J., Luhman, K.L., Muzerolle, J. & Young, E. 2007, AJ, 134, 411 
Muench, A. A., et al. 2003, AJ, 125, 2029 
O'Dell, C.R. 2001, ARAA, 39, 99 

Ouellette, N., Desch, S. J., Hester, J. J. & Leshin, L. A. 2005, in Chondrites and the 
Protoplanetary Disk, ed A. N. Krot, E. R. D. Scott & B. Reipurth (San Francisco: 
ASP Conference Series), 527 

Ouellette, N., Desch, S. J. & Hester, J. J., ApJ, 662, 1268 

Palla, F. & Stahler, S. W. 2000, ApJ, 540, 255 

Palla, F. & Stahler, S. W. 2001, ApJ, 553, 299 

Palla, F., Randich, S., Pavlenko, Y. V., Flaccomio, E. & Pallavicini, R. 2007, ApJ, 659, L41 
Preibisch, T. & Zinnecker, H. 1999, AJ, 117, 2381 

Quitte, G., Halhday, A. N., Markowski, A., Meyer, B. S., Latkoczy, C. & Giinther, D. 2007, 
ApJ, 655, 678 

Quitte, G., Bizzarro, M., Markowski, A. & Latkoczy, C. 2006, MAPS, 41 (Supp), #5236 



Rauscher, T., Heger, A., Hoffman, R. D. & Woosley, S. E. 2002, ApJ, 576, 323 

Reacli, W. T., et al. 2004, ApJS, 154, 385 

Regelous, M., Elliott, T., Coath, CD., 2007, GCA, 71 A825 

Richling, S. & Yorke, H.W. 2000, ApJ, 539, 258 

Rho, J., Reach, W. T., Lefloch, B., & Fazio, G. G. 2006, ApJ, 643, 965 

Romano, D., Chiappini, C., Matteucci, F. & Tosi, M. 2005, A&A, 430, 491 

Russell, S. S., Hartmann, L., Cuzzi, J., Krot, A. N., Gounclle, M. & Wcidcnschilling, S. J. 
2006, in MESS2, ed D. S. Lauretta & H. Y. McSween (Tucson: AUP), 233 

Russell, S. S., Gounelle, M. & Hutchison, R. 2001, Phil. Trans. Royal Society A, 359, 1991 

Salpeter, E.E. 1955, ApJ, 121, 161 

Sandell, G. & Knee, L. B. G. 2001, ApJL, 546, L49 

Scalo, J. M. 1986, Fundamentals of Cosmic Physics, 11, 1 

Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 1992, A&AS, 96, 269 

Shukolyukov, A. & Lugmair, G. W. 1993, EPSL, 119, 159 

Smith, N. & Brooks, K. J. 2007, MNRAS, 379, 1279 

Smith, N. 2006, MNRAS, 367, 763 

Smith, N. et al. 2000, ApJ, 532, L145 

Smith, N., Bally, J. & Morse, J. A. 2003, ApJ, L105 

Srinisavan, G., Goswami, J. N. & Bhandari, N. 1999, Science, 284, 1348 

Storzer, H., & HoUenbach, D. 1999, ApJ, 515, 669 

Sugerman, B. E. K., et al. 2006, Science, 313, 196 

Sugiura, N., Miyazaki, A. & Yin, Q. Z. 2006, Earth Planets Space, 58, 1079 
Throop, H. B., & Bally, J. 2005, ApJ, 623, L149 

Tachibana, S., Huss, G. R., Kita, N. T., Shimoda, H. & Morishita, Y. 2006, APJ, 639, L87 



-28- 



Tachibana, S., Huss, G. R. & Nagashima, K. 2007, LPSC, 38, 1709 
Tassis, K., & Mouschovias, T. C. 2004, ApJ, 616, 283 
Urey, H. C. 1955, Proc. Nat. Acad. Sci. U.S., 41, 127 

Wadhwa, M., Amelin, Y., Davis, A. M., Lugmair, G. W., Meyer, B. S., Gounelle, M. & 
Desch, S. J. 2007, in Protostars and Planets V, ed B. Reipurth, D. Jewitt & K. Keil 
(Tucson: University of Arizona Press), 835 

Wang, W. et al. 2007, A&A, 469, 1005 

Wasserburg, G. J., Busso, M., Gallino, R. & NoUett, K. M. 2006, Nuclear Physics A, 777, 5 
Weaver, R. et al. 1977, ApJ, 218, 377 
Williams, J.P. & Gaidos, E. 2007, ApJ, 663, L33 
Woosley, S. E. & Weaver, T. A. 1995, ApJS, 101, 181 
Wren, Z., O'Dell, C.R. 1995, ApJ, 438, 784 



This preprint was prepared with the A AS IM^jX macros v5.2. 



- 29 - 




Fig. 1. — Evolutionary timescales of massive stars. Schaller et al. (1992) fitted tfieir models 
with the expression log(tsN) = l-4/(logMsN)^'^ where tsN and Msn are respectively the 
evolutionary timescale (in Myr) and the mass of the considered star (in Mq). More recent 
work by Romano et al. (2005) confirm the findings of Schaller et al. (1992). 
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Fig. 2. — Sketch illustrating the fundamental problems met by the idea of short-lived ra- 
dionuclides injection in a protoplanetary disk by a nearby supernova. In this cartoon, the 
evolution of a massive star (large circle) and several low-mass stars within molecular gas 
(light grey) are shown. In panel (a) is depicted the onset of star formation assumed to be 
coeval for all stars. During that stage, only dense cores are present. Note that if there is 
any spread in the onset of star formation, it would correspond to a delayed formation of the 
massive star (see text). One million year later (panel (b)), the massive star has carved a 
~ 0.3 pc wide HII cavity. Disks which are within this cavity are evaporating as depicted 
by their cometary appearance similar to that observed in Orion protoplanetary disks. Star 
formation is less vigorous than at T = Myr. At T = 2 Myr (panel (c)), the HII cavity has 
expanded to ~ 1 pc (see text). The disks which are inside the cavity are totally photoevapo- 
rated (Johnstone, Hollenbach & Bally 1998). Low-mass stars that formed at T = Myr now 
have planets. Very few new stars are forming (see text and Figure H]). When the massive 
star goes supernova (panel (d)), the stars in the uttermost vicinity of the supernova do not 
have disks, while disks further away have already dissipated to form planets. At this stage, 
star formation is virtually non-existent because molecular gas has largely been cleaned by 
the massive star wind and UV emission. A left-over, evanescent, parcel of molecular gas is 
shown in light grey. Injection of SRs by the supernova in a late-formed disk (marked with 
an arrow) is an extremely improbable event. An indicative scale would be given by the size 
(few tenths of parsec) of the HII cavity in panel (b). 
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Fig. 3. — Probability that the most massive star in a cluster containing A^^, stars has a mass 
-^SN = 40 M0 and Msn = 120 Mq (equation (jl])). Note that the maximum probability for 
a 40 M0 star to be the most massive star occurs for a cluster of ~ 4000 stars, slightly larger 
than the ONC. 
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Fig. 4. — Star-formation rate versus age for a suite of well studied molecular clouds. Cham 
and p Oph stand for Chamaeleon I and p Ophiuclii respectively. The thick dashed line 
is the fit of the data which minimizes the exponential decrease of the star formation rate 
(see section 14.41) . It gives a decay time constant of r = 0.56 Myr. The age of a given 
star-forming region depends on the mass of the stars and their position on the Hertzsprung- 
Russell diagram, among other parameters (Palla & Stahler 2000; Hartmann 2001). For p 
Ophiuchi, we take the commonly accepted age of ~ 1 Myr (Preibich & Zinnecker 1999; Palla 
& Stahler 2000). For other molecular clouds, we adopt the age compilation as well as their 
associated error bars given by Haisch, Lada & Lada (2001). These ages are compatible with 
other, more recent, estimates in the literature (e.g. Luhman 2004; Muench et al. 2007). The 
errors take into account the error in the mean of the source age derived from a given set 
of Pre Main Sequence (PMS) tracks (Haisch, Lada & Lada 2001). Errors due to the use of 
different PMS models are of the order of 1 Myr and could result in a shift of the data without 
changing the relative age sequence (Haisch, Lada & Lada 2001). Star formation rates for 
individual molecular clouds are given by Preibisch & Zinnecker (1999) without associated 
error bars. Note that star formation rates are not available for older clusters because there 
are no observed star-forming molecular clouds older than a few Myr (Elmegreen 2000; Lada 
& Lada 2003; Hartmann, Ballesteros-Paredes & Bergin 2001). 
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Fig. 5. — Fraction of low-mass stars younger than 1 Myr (/yso) present in the cluster as a 
function of time t (equation [10]). /yso is shown for the 3 expressions considered for the star 
forming rate, ip{t), and discussed in the text. 
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Fig. 6. — Probability of protoplanctary disk enrichment in SRs at the level observed in the 
solar system (Psr(-^*)) as a function of the cluster size (A^*), and calculated for the different 
expressions of the star forming rate ip{t). The curves were obtained using the ONC stellar 
surface density of Sqnc = 211 stars/pc^. A similar curve obtained by Williams & Gaidos 
(2007) is shown for comparison (with S3 = 211 stars/pc^). 
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Fig. 7. — Probability for any protoplanetary disk attached to a young low-mass star to 
receive SRs from a nearby SN at the level observed in the solar system. The vertical line 
denotes the surface density of the ONC, while the grey circle represents our preferred model, 
with the Orion surface density and a linear decrease of ip{t). 
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Fig. 8. — Number of stars with mass larger than 40 Mq as a function of cluster size (calculated 
from equation (^). Multiple supernovae 2) capable to inject SRs in a cluster occur only 
in clusters containing more than 10'^ stars. 



